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Abstract 

Using the method of QCD light-cone sum rules, we calculate the B — * tttt hadronic 
matrix elements with annihilation topology. We obtain a finite result, including the 
related strong phase. Numerically, the annihilation effects in B — > tttt turn out to be 
small with respect to the factorizable emission mechanism. Our predictions, together 
with the earlier sum rule estimates of emission and penguin contributions, are used 
for the phenomenological analysis of B — > tttt channels. We predict a A 7 = 1/2 
transition amplitude which significantly differs from this amplitude extracted from 
the current data. 



1 Introduction 



Charmless B decays play an increasing role as a possible window to physics beyond Stan- 
dard Model (for a review see. e.g., [1]). In particular, the decays of the type B — > nn and 
B — > Kn,KK are actively investigated using various approaches [2-12]. 

One promising method to evaluate charmless nonleptionic B decays is QCD factor- 
ization (QCDF) [2,3], which to leading order in l/m b and a s (m b ) yields the well known 
formulae of naive factorization. However, even when QCD corrections and estimates for 
the penguin contributions are included, QCDF is only marginally consistent with the data 
for B — > 7T7T. As a more general analysis of data shows, one indeed needs additional 
contributions to the decay amplitudes which correspond to the AI = 1/2 piece of the 
Hamiltonian. Comparing with different topologies of the decay diagrams, this means an 
enhancement of either the penguin or annihilation contributions. Especially interesting is 
the weak annihilation (or weak exchange) of the b quark and the light antiquark in the 
B meson. In QCDF, formally suppressed as l/mb, the annihilation contribution is diver- 
gent. Fits of QCDF to data [4], with the annihilation part replaced by finite parameters, 
yield large effects. Importantly, the fits themselves cannot clearly distinguish between the 
annihilation and penguin mechanisms. In order to assess the relative importance of the 
annihilation contribution, one needs a different approach in QCD which allows to calculate 
B — > 7T7T hadronic matrix elements with various topologies. 

In the present paper we will employ the technique suggested in [13]. It is based on 
light-cone sum rules (LCSR) [14], the method which adapts the general idea of QCD sum 
rules [15] for the amplitudes of exclusive hadronic processes. One well known application 
of LCSR is the calculation of the B — > 7r form factor [16]. In the sum rule approach, as 
opposed to QCDF or to the perturbative hard-scattering QCD approach (PQCD) [6], one 
does not directly represent the hadronic matrix element in terms of quark-gluon diagrams. 
The QCD calculation takes place for the correlation function, a more general object, where 
the separation of short and long-distance parts into hard-scattering amplitudes and pion 
distribution amplitudes (DA's) respectively, is possible "by construction", that is, due to 
an appropriate choice of kinematical variables. Both hard and soft gluon effects are sys- 
tematically included in this calculation, contributing to different terms of the light-cone 
operator-product expansion (OPE). The hard gluon exchanges enter the hard-scattering 
amplitudes, whereas the soft gluon effects are represented by quark-antiquark-gluon DA's 
of the pion. The hadronic matrix element appears as a part of the hadronic dispersion 
relation for the correlation function. Matching this relation to the result of the QCD 
calculation, one employs quark-hadron duality to separate the ground-state B — > nn am- 
plitude from the background of excited hadrons. The advantage of this technique, as will 
be explained in more details further, is a possibility to associate different decay topologies 
in the hadronic matrix element with the corresponding diagrams in the OPE, hence the 
relative contributions of various operators and topologies to the decay amplitude can be 
estimated. Importantly, the calculation of B — > nn amplitudes takes place in full QCD at 
finite m h and one uses the same input as in the sum rule for the B — > n form factor. The 
latter provides also the factorizable part of the B — > nn amplitude. 
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Penguin topology contributions to B — > 7T7t have been estimated using LCSR already in 
previous papers [17,18], while annihilation topologies are notoriously difficult to calculate. 
This is due to the fact that the factorizable contribution of the current-current operators 
Oi ; 2 with annihilation topology vanishes. Hence, in a sum rule calculation, the leading 
nonfactorizable annihilation with one hard-gluon exchange already corresponds to a set of 
two- loop diagrams with several scales (b quark mass and external momenta). 

In the present paper we employ a simplified method for the annihilation via hard 
gluons, that is, both final state pions are replaced by DA's and only the initial state B 
meson is interpolated by an appropriate current. Importantly, the propagators appearing 
in the resulting Feynman diagrams carry sufficient virtuality to be treated perturbatively. 
Although this is a modification of the original idea presented in [13], it is a useful approach, 
since it yields a result free of infrared divergences. This result may be compared to what 
is obtained from QCDF, where an infrared regulator needs to be introduced, reducing 
basically the whole annihilation contribution to a nonperturbative parameter of order 1/mj,. 
On the other hand, the annihilation with soft-gluon exchange which is not accessible in 
QCDF, is calculated here within the standard procedure [13] leading also to a finite answer. 
After adding the hard and soft-gluon contributions from LCSR, we obtain the main result: 
the finite matrix element of the operator 0\ with annihilation topology, including its phase. 
We predict this effect to be numerically small. 

Furthermore, we investigate the annihilation contributions via quark-penguin operators. 
While 03,4 have the same V — A structure as 0± t 2, so that the factorizable annihilation 
vanishes, the operators 5j6 contribute through two different types of contractions. One 
of them has a V + A content and also vanishes in the factorizable annihilation, whereas 
the other one, with a S ± P structure allows for a factorizable B — ► mi transition with 
annihilation topology. This contribution reduces to a separate nonperturbative object, the 
pion scalar form factor at timelike momentum transfer m 2 B . In QCDF and PQCD, the 
factorizable annihilation was taken into account only with a perturbative gluon exchange 
between the final state quarks, corresponding to the pion form factor in 0(a s ). The method 
of LCSR allows to obtain the zeroth order in a s , that is the "soft" (end-point) part of the 
scalar pion form factor. We calculate this part, with the same approach as for the e.m. 
(vector) pion form factor in [19,20]. The resulting hadronic matrix element is large, due 
to the chirally enhanced factor. Still this effect alone cannot produce a large AI = 1/2 
amplitude in B — > -kit, because of the small Wilson coefficients. 

Finally, we evaluate the B — > tt7t decay amplitudes using LCSR predictions and includ- 
ing all calculated nonfactorizable effects with the emission, penguin and annihilation topol- 
ogy. Since these effects are generally small, the discrepancy between the B° — > 7r + 7r~, 7r°7r° 
observables calculated in the factorization limit and the current experimental data remains. 
One may encounter a situation similar to K — > titt, with its long-standing problem of 
AI = 1/2 rule. 

The paper is organized as follows: In section |2] we summarize the current status of 
phenomenology of B — > inr amplitudes, starting from the isospin decomposition and naive 
factorization and comparing them with the data. In section El we discuss nonfactorizable 
effects, representing the B — > tttt decay amplitudes in terms of hadronic matrix elements 
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of effective operators with different topologies. In section HI we derive the LCSR for the 
hadronic matrix element of an effective operator with a given topology. In sections 03 
and El we present our new results for B — > tttt annihilation with perturbative (hard) and 
nonperturbative (soft) gluons, respectively. In section we obtain the sum rule for the 
pion scalar form factor which determines the specific factorizable annihilation contribution 
of the operators 0§q. In section |H1 we perform the numerical analysis and present the 
LCSR prediction for the annihilation contributions to B° — ► tt + tt~ . Furthermore, we add 
all calculated contributions and present our numerical predictions for the branching ratios 
and direct CP-asymmetries in all three B — > tttt channels. In section E3 we analyze our 
result in the limit m& — > oo and comment on the annihilation mechanism in QCDF and 
PQCD. We conclude in section [TUl The appendices contain some expressions used in the 
paper. 



2 Phenomenology of B — > tttt amplitudes 

Throughout this paper we adopt isospin symmetry. We also neglect the effects of the 
electroweak penguin operators, so that the effective weak Hamiltonian for B — > tttt has the 
following expression: 

H eS = ^{A u (ci 01 + c 2 O u 2 ) + A c (d 0\ + c 2 2 c ) 

6 

+ (A u + A c )[^c i O i + c 8 p0 8g ]} + h.c. . (1) 

i=3 

The CKM factors are defined as A p = V P bV* d (p = u,c,t), and we use the CKM unitarity 
replacing — X t by \ u + A c . Hereafter we suppress for brevity the scale dependence in the 
Wilson coefficients q, which is supposed to be compensated by the scale-dependence of 
the hadronic matrix elements of the effective operators Oj. The current-current operators 
entering Eq. (JJJ are 

of = (<Ii» (pm) = l -oi + idi , o\ = (pr, P ) (dr»b) = l -oi + 2d{, (2) 

where p = u,c, = 7 M (1 - 75) and 



o\ = (dr^p)(pr^b) , 01 = (pr^p)(dr^b) (3) 



with Tr(A a Aft) = 28^. The color Fierz transformation allows to use, instead of the combi- 
nation C1OI + C2O2, either (ci + c 2 /3)Of + 2c 2 0f or the opposite one with 1 <->• 2. In leading 
order the operators with color-neutral currents factorize and nonfactorizable contributions 
start from two-gluon exchanges, which we will neglect. The color-octet currents yield non- 
factorizable effects starting at a one-gluon level. These effects will be systematically taken 
into account. 
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To complete the definition of H cS , we specify the quark penguin operators: 

Os = X)(/T/J)(<fln) = \0 A + 20 4 , 4 = J2(dTJ)(fT%) = io 3 + 20 3 , 

/ / 

5 = Y^Hi 1 + 75)/) (dT"&) = ^ 6 + 20 6 , 
/ d 

6 = -2 £(d(l + 75)/) (/(I - 7fe)6) = ^ 5 + 20 5 , (4) 

where / — u,d, s, c, 6 and we again use the color Fierz decompositions introducing the 
operators with color-octet currents 03,4,5,6 obtained from 03,4,5,6 , respectively. Finally, the 
chromomagnetic quark-gluon penguin operator is: 

°s 9 = -^m b da, v (l + 75 )G^6 . (5) 

Turning to the phenomenology of B — > nn we begin with quoting the results of the 
current measurements [25]: 

BR(B + -> tt+tt ) = (5.5±0.6)xl0 -6 
BR(B° ^ 7T + 7i-) = (5.0±0.4)xl(T 6 

BR(B° -> ttV) = (1.45 ± 0.29) x 1(T 6 , (6) 

where only the first branching ratio is compatible with the expectations of factorization. 

The problem may be analyzed in terms of the usual isospin decomposition. The effective 
Hamiltonian of the Standard Model consists of two parts with AI =1/2 and AI = 3/2, 
resulting in two reduced isospin amplitudes A and A 2 in B — > nn, which correspond to 
the pions in the 1 = and 1 = 2 final states, respectively. One obtains the following 
decomposition for the amplitudes 1 

A (B- - tt-tt ) = (tt-tt !^^-) =A A 2 , 

-> 7T + 7T-) = (TT + TT-liJeffl.B ) = A 2 + A , ( 7 ) 

- ttV) = (TrV^efflfi ) =2A 2 -A , 
from which the well known isospin relation [21] is obtained: 

A(B° -> 7T 7T ) = V^A(S" -> 7T-7T ) - A(B° -> 71 + . (8) 

In the above we use the same convention for the amplitudes, as in Ref. [2,3], including the 
statistical factor 1/2 in the branching ratio for B° — > -r°-k°. 

1 Throughout this paper we consider, for definiteness, B~ and B° decay amplitudes, whereas all quoted 
branching ratios for B + and B° are, as usual, CP-averaged. 
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From Eq. (JJJ) one obtains the ratio of the moduli of Aq/A 2 in terms of the decay rates 



\Ao\ 



BR(B° -> 7T°7i°) + BR(B° -> 7T+7T-)] r B - 



2. (9) 



-> 7T-7T ) 

Using Eq. © (neglecting CP asymmetries) and t b -/t b o = 1.075 ± 0.009 [25], one gets 

\A \ 



\Ao 



1.33 ±0.31. (10) 



Employing the tree-level emission graphs and retaining only the current- current opera- 
tors (the naive factorization limit) we may obtain a first insight into the anatomy of these 
decays. In fact, this simplified way leads practically to the same conclusions as the full 
calculation in the framework of QCDF. First of all, the decay B~ — ► 7r~7r° is well described 
in the factorization limit, where 

V2A(B- -> tt-tt ) = X u ~ (ci(/i) + C 2 (/i)) A nw (11) 

with the usual notations for the Wilson coefficients ci )2 and the factorizable B — ► ttti 
amplitude 

^ = »|/Xk)K-^). (12) 

In the above, f n = 131 MeV is the pion decay constant and f Bn (q 2 ) is the scalar B — > 7r 
form factor. Hereafter, we neglect in the amplitudes, retaining it only in the ratio 
fi* = ml/(m u + m d ). 

In the factorization approximation, it is sufficient to use the leading-order (LO) values 
of the Wilson coefficients. We vary their renormalization scale within m^/2 < /i < mj; for 
illustrative purpose we also put /x = Mw- Using the LCSR prediction [16] for the B — » n 
form factor: 

- /L(0) = /^(0) = 0.26 ± 0.05 , (13) 

which is explained below in sect. |HJ and taking \V u b\ = (4.22 ± 0.11 ± 0.24) x 10 -3 from 
Ref. [22] (adding the errors in quadrature) we obtain from Eq. (jllj) 

{{5.7±H ± 0.7) x 10~ 6 , /i = m b /2 (a = 1.169, c 2 = -0.361) 
{QAtli ± °- 8 ) x 10-6 , ^ = wit (ci = 1.108, c 2 = -0.249) 
(8.711-2 ± 1-1) * 10~ 6 , n = M w ( Cl = 1, c 2 = 0) , 

(14) 

where the errors reflect the uncertainties of the form factor and of \V u b\, respectively. The 
scale-dependence is mild, as expected. 

Another way to check the validity of the factorization approximation for this channel 
(independent of \V u b\ and the value of f B7V (0)) is provided by the ratio of the B + — ► 7r + 7r° 



5 



and B° — > tt l + V[ widths. In factorization approximation: 

BR(B+ tt+tt ) 27r 2 |Kd| 2 (ci(/i) + c 2 (/i)) f 2 m% , Tb+ . 



(15) 



where E n = (m 2 B + m 2 — q 2 )/(2m B ) and N B7T {q 2 ) = /^(^V/st^O) is the shape of the 
form factor. The recent measurement [23] of the B — > nlu decay distribution, fitted to 
the parameterization [24], N Bw (q 2 ) = [(1 — q 2 /m B *)(l — a BlT q 2 /m 2 ^)}' 1 , yields olb-k — 
0.61 ± 0.09. Using this value, and the average experimental number BR(B° — > tt^Z+z/^) = 
(1.36 ± 0.11) x 10~ 4 [25], we obtain at n = m b /2 [m b ; M w \: 

BR(B + -> n + n°) fact = (3.6 ± 0.3 ± 0.4) x 10" 6 

[(4.1 ± 0.3 ± 0.5) x 10~ 6 ; (5.6 ± 0.5 ± 0.6) x 10~ 6 ] , (16) 

where the first error originates from the semileptonic branching ratio, and the second one 
from the slope parameter a Bn . The result is again in the ballpark of the experimental 
interval (JBJ). Hence we conclude that the amplitude A 2 may be estimated correctly by 
naive factorization. 

Furthermore, using naive factorization, we may express the ratio A /A 2 in terms of the 
Wilson coefficients, since all hadronic matrix elements will drop out. Using Eq. (fTT|) and 
the analogous relation 

A(B° - 7T + 7T-) = K(cM + (17) 



and comparing with the decomposition ((7j) we obtain 



Aq _ 5 / c x (/i) - c 2 (//)/5 \ 
A 2 ~ Ay ci(/x) + c 2 (//) J 



1.92,/i = m b /2 
1.68 , [i = nib 
1.25, n = M w . 



In fact, this expression depends quite strongly on the scale, showing that naive factorization 
for B° modes is compatible with the data only for a large scale of 0(niw), which seems 
unrealistic. 

A more distinct disagreement is revealed between the ratios of B° and B + partial widths 
calculated in the naive factorization at fi = nib/ 2 [mb] Mw}'- 



BR(B° -> tt+tt-) BR(B 



_ ^0 



1.77 [1.49; 1.05], - D / p ^ ~ 10~ J [0.010; 0.058] , (19) 



BR(B+ TT + TT ) ' £#(£ + -> TT + TT 

and the same ratios obtained from the experimental results (JHJ). In fact, the calculated B — > 
TT Tr° width is too small even at /i = m^. We conclude that the naive factorization picture 
misses an important part of the amplitude Aq which interferes destructively (constructively) 
with A 2 in A(B° — > tt+tt - ) (A(B° — > tt tt )). If there were large nonfactorizable gluon 
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corrections to the emission topology, they would have influenced both A and A 2} violating 
the abovementioned agreement for the B~ — > 7t~tt° channel. Hence, the missing isospin 
zero amplitude should be searched for within the contributions of non-emission topologies 
for the current-current operators and/or in the contributions of the penguin operators, an 
opinion shared by many recent analyses of these decays. 

3 Beyond Factorization 

In Fig. ^ we schematically represent different quark topologies contributing to the am- 
plitude of B° — > 7r + 7r~, the channel of our interest. Under topology we understand the 
way to contract the valence quarks (antiquarks) of the initial B and final mesons with the 
antiquarks (quarks) from the effective weak Hamiltonian. This concept is meaningful for 
all those methods where the valence quark structure of mesons is well defined, either via 
the meson distribution amplitudes (DA's) as, e.g. in QCDF, or via both DA's and inter- 
polating currents, as in the correlation function for LCSR. Denoting the valence spectator 
antiquark in B by q s and the quarks emitted in the b-quark decay by q\,q2,q3 we define 
emission as the part of the decay amplitude where all four quarks and antiquarks end up 
as the valence quarks of the final mesons (Fig. ^i). Correspondingly penguin is the part 
where q s and only one of qi,q2 belong to the final mesons (Figs. E),c). The remaining 
two possibilities are annihilation (Fig. [Hi ) or penguin annihilation (Fig. Hk), where either 
<7~3 and the one of qi , q<i or none of the quarks are among the valence quarks of the final 
mesons, respectively. 2 

In the isospin symmetry limit it is sufficient to investigate in detail the B° — > 7t + tt~ 
amplitude. It has a single 1 = 2, 0(X U ) part with emission topology, common with 
the B~ — > 7T~7r° amplitude, and in addition contains many different 1 = contributions 
proportional to both X u and A c . The B° — > 7r°7r° amplitude is then obtained by simply 
using the isospin relation (jHJ). 

The complete isospin decomposition of the B° — > 7t + tt~ amplitude following from 
Eq. (0) can be cast in the following form: 

A{ B° -> 7T+7T-) = X u (4^' 2) + 4^ 2) ) + X C A^' 2) + (X u + X C )A^ 3) , (20) 

where the upper indices in Aj indicate the contributing operators and the lower index the 
isospin. Correspondingly, 

V2A(B- -> tt-tt ) = A u (34"' 1,2) )- (21) 

Each separate amplitude A® in Eqs. (j20l21j) contains a sum over hadronic matrix ele- 
ments with different topologies (T): emission (E), penguin (P q ,P c ,P b ), annihilation (A) 

2 Note that drawing generic quark-line diagrams, one can always start from the emission topology 
and then, merging the quark- antiquark (spectator) lines in the final state, end up with an annihilation 
mechanism; according to our classification, this mechanism still belongs to annihilation. 
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(a) emission 



(b) penguin 



(c) charming penguin 




(d) annihilation 




(e) penguin annihilation 



Figure 1: Different quark topologies in B° 
lines the W-boson. 



7r + 7r ; double lines denote the b quark, wavy 



and penguin annihilation (P q A,P c A,PbA) , where q = u,d,s. According to Ref. [26], this 
subdivision can be systematically done in a scheme- and scale-independent way. It is clear, 
that the number of independent matrix elements of the type (7r + 7r~|Oj|i? )T is less than 
the number of operators because certain penguin operators have the same quark struc- 
ture as the current-current operators. For the penguin topologies we will also neglect the 
differences between quark loops with q = u,d, s. 

Let us first discuss the 1 = 2 part in Eqs. (ffilj) and (j2H), which is relatively simple. 
Taking into account the nonfactorizable emission correction, where only 0\ contributes in 
the one-gluon approximation, we introduce the ratio: 



^("7771") 



(tt + tt 


0\ 


B°) E 




0\\B*) E 



(22) 



Note that in the adopted approximation the matrix element standing in the denominator 
is factorizable: Gf/V2(tt + 7i^\O'i\B ) e = A n7T . We obtain: 



.4 



0,1,2) 



^[ci + c 2 ) + -(ci + c 2 )r E 



(tttv) 



A 



(23) 



Accordingly, the 1 = part generated by the operators 0\ 2 has the following decomposi- 
tion: 



.4 



(m,1,2) 



5 ISp. 



C 2 



-[ Cl -2c 2 )r E 



{nw) + 2c 1 (r { ™ ) + r { ™ ] 



A 



(24) 
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where the relative contributions for penguin and annihilation topologies 3 similar to Eq. ()22|) 
are defined as follows: 



(tttt) = ^TT'TT \U5\£r)p u frw) _ (7T+7T \0%\B°) A 

Tp « " (n+n-\Of\B°) E ' Ta (n+n-\0^\B°) E ' 1 ' 

Hereafter, we will use a more generic notation q = u,d,s } so that rp^' = r p ^ in 

accordance with our approximation for the light-quark penguin loops. Furthermore, the 
0\ 2 operators with c-quark contribute to Eq. (f271jl only in the penguin topology ( "charming 
penguins"): 

.(0,1,2) Q (tttt) , (tttt) {^ + ^\Ol\B Q )p c 

A\ = 2c\r p ' where r\> ' = - — ■ — -=—. — . (26) 

The remaining piece of the decomposition (f2"U|) containing the hadronic matrix elements 
of quark-penguin operators with various topologies, normalized to the factorizable part, is 
more complicated: 



(>3) 



C3 . / (tttt) (tttt) (tttt) (wtt) \ . n ( n -(nn) -(tttt) -(tttt) (tttt) 

c 4 + — + 2c 3 {r E +r Pq + r Pb + r\ >) + 2c 4 \3r), + r p + r\ + 2r v 



2/^7T ( C 5 \ n (7T7T,6) . r, / -(7T7r) -(tVIt) -(tV1t) (7T7r,5) 

( c 6 + v ) + 2c 5 r E + 2c 6 1 3r> +r Pc J + r Pb ' + 2r\ 



1 / 1 r>(7r7r,6) . (7T7t,6) . eff (tttt) 

V 6 37 A +2c5rA +% r s 9 



Some of the -parameters entering the above equation have already been defined: 
we use the fact that certain quark-penguin and current-current operators coincide. The 
parameter r p determines the relative contribution of the 6-quark penguin topology. It 

is defined as in Eq. (126(1 with c — ► b. Furthermore, the notation f P b is introduced to 

distinguish the penguin contractions of the operators 4j6 (0 3j5 ) from from those of 0i 3 
(^2,4). In the NDR scheme used here, as in [2,3], the quark loop factors for the two 
contractions differ by a constant. The loop factor entering LCSR for is given in 



Eq. (4) of [18]. To obtain the corresponding factor for r^ 71 "^ one simply has to subtract 
1/6 from this expression. A few terms in Eq. (|27|) are generated by the effective operator 
Oq with (S + P) ® (S — P) structure. First, the factorizable emission contribution of 
the w-quark part of this operator, Og = — 2(d(l + ys)u)(u(l — 75)6), acquires a "chirally- 
enhanced" factor fi^/nib, whereas its nonfactorizable part is described by an additional 
parameter: 

r (^,6)_ _ _ (2S) 



(7T + 7T 




B°) E 


(7T + 7T - 


Ot\B») E 



In this paper, we neglect the penguin annihilation (PA) topologies. They can be simply added to the 

(tttt) 
PA 



general decomposition by introducing the corresponding Tp^) ratios, but from LCSR we expect them to 



be small; this mechanism is also neglected in QCDF. 
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Second, there is a specific factorizable annihilation contribution due to the <i-quark part of 
the same operator Oq = — 2(d(l + j 5 )d)(d(l — 7s)&) which is expressed in terms of the B 
meson decay constant {Q\d^^b\B) = —im^fB/mb multiplied by the pion scalar form factor: 

(vr + (pi)7r-(p 2 ) | dd | 0) = F*({ Pl +p 2 f) , (29) 

corresponding to the transition of the scalar and isoscalar quark-antiquark current into a 
two-pion state with the invariant mass squared (P1+P2) 2 = ^b- The parameter in Eq. (|2Tj) 
determining the annihilation via Og is factorized as: 

R (^6) (n+n-\Ol\B°} A 2f B Ff(m%) 



The nonfactorizable annihilation correction of the color-octet counterpart of this operator 
in Eq. (}2Tj) is parameterized as 

In addition to this S ± P contractions, the nonfactorizable annihilation contribution from 
the V+A contraction is defined in analogy to Eq. (j31j) with (6 — > 5). Finally, the parameter 
in Eq. (|2*7jl describing the contribution of the gluonic penguin operator (with the penguin 
topology) is: 

,("■"") _ \" " ! ~ v / <" ' 1 (32) 



(/T + 7r 


Os g \ 


B°)p g 


(7T + 7T - 


\OZ\B°)e 



In the NDR scheme it should be multiplied by cJJ = c Sg — c 5 (see 

For convenience, in Appendix A we present the relations of the parameters intro- 
duced above to the effective coefficients ai and 6, used in QCDF [2,3] to encode the non- 
factorizable effects in B —>■ tttt. Eq.(|2()jl can also be converted into a typical decomposition 
in terms of "tree", "color-suppressed" and "penguin" amplitudes used in the CP-analysis 
of charmless decays, where the separation to X u = |A u |e~ i7 and A c parts is made explicit: 

A(B° -> 7T+7T-) = e-^T nw + P m , y/2A(B~ -> t^tt ) = e"^(T_ + C^) . (33) 

Comparing the above with Eq. (f2*U|) one reads off: 

T — I \ if 1 4(^1.2) , a(>3)\ 

1 7T7T — I A U I I ^2 *T" ^0 ' ^0 ) 1 

PlTTT = A C (4 ' ' ^ + ' Pirn "I" C*7r7r = | A u | ^3A^ ^ . (34) 

To analyze the B — > titt amplitudes in terms of separate isospin contributions, one 
needs the numerical values of the hadronic parameters entering Eqs. (|23 |) . ()24jl . (|2fi|) and 
(J2*7J) . Some of them have already been estimated using LCSR. For the nonfactorizable 
emission entering through we will partly use the QCDF result. The calculation of 
the unknown annihilation parameters and R^' 6 ^ from LCSR is the main issue of this 
paper. 
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4 Derivation of LCSR 



The method we apply is basically the one developed in [13], (see also [17,18]) with some 
important modifications which will be explained in this section. To demonstrate the deriva- 
tion of LCSR for a generic hadronic matrix element (7r + 7r~|0|-B°)T, let us choose for def- 
initeness the combination of current-current operators O = c\0\ + C2O2 at a fixed scale, 
considering it as a superposition of two local operators. 
One starts from the correlation function 

Fi°\p,q,k) = -Jd 4 x e-^'^J d 4 y e ^"^(0|T {j${y)0*(Q)ji B \x)} \*-(q)), (35) 

where = u r ) a r )^ ) d and j'5 = mjn^d are the quark currents interpolating the pion and 
the B meson, respectively. The momentum k is artificial and will vanish in the final sum 
rule. It is introduced in order to have two independent kinematical variables in the B and 
7T7T channels. 

The correlator (J35|) can be decomposed into four different Lorentz structures, 
Fi° ] = (p - fc) a F<°> + q a F[ 0) + k a F^ + e^gVFf) 

of which we use only the first. Concerning the kinematical variables, we put q 2 = m 2 = 
and choose p 2 = k 2 = for simplicity. The remaining invariants are (p — k) 2 , (p — q) 2 and 
P 2 = (p — k — q) 2 . In the domain where all three variables are spacelike and large, all 
distances are close to the light-cone, x 2 ~ y 2 ~ (x — y) 2 ~ 0, and the correlation function 
can be calculated by perturbatively expanding the T-product of operators. In this way, the 
correlation function is expressed in a usual form of hard-scattering amplitudes convoluted 
with pion DA's of growing twist. 

For a given operator O in the correlation function (|35jl . various contractions of the quark 
fields are possible, leading to diagrams with different topologies. Collecting the lowest-order 
contributions to OPE, we easily recognize diagrams with the emission (Fig. EJ), penguin 
(Fig. I2J), annihilation (Fig. HJ) and penguin- annihilation (Fig. EJ) topologies. The sum of all 
diagrams calculated at a definite order in OPE, will be matched to the dispersion relation 
where the ground-state contribution contains the hadronic matrix element (7r + 7r~|0|i? ). 
If one retains only diagrams with a topology T in the OPE, the sum rule result (within 
adopted accuracy) can be interpreted as (7r + 7r~|O|.B )T- 4 This was actually done for the 
emission topology in [13] where the diagrams in Fig.[2]have been investigated. It was shown 
that retaining only the diagram of Fig. Eh, without gluons connecting the light-quark loop 
and the heavy-light part, one reproduces the result of naive factorization. Importantly, the 
gluons which do not violate factorization in this diagram can be added arbitrarily. Alto- 
gether, one obtains the product of the LCSR for the B — > 7r form factor and the two-point 

4 Since we are studying only the leading-order effects here and use a fixed scale, the scale- and scheme 
dependence as well as the mixing effects between separate operators in -ff e ff, remain beyond our scope. To 
account for these effects one has to consider scale- and scheme invariant combinations of matrix elements 
with different topologies as explained in [26]. 
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(c) 

Figure 2: Diagrams corresponding to the emission topology in the OPE of the correlation 
function \35)) : (a) factorizable; (b) with nonfactorizable hard gluon (six diagrams) ;(c) non- 
factorizable soft gluon (two diagrams). The solid, double, dashed, wavy lines and the square 
denote the light quarks, b quark, gluon, external currents and the weak vertex, respectively. 
The shaded ovals denote the pion DA 's. The crosses indicate how gluon lines are attached 
in the other possible diagrams. 




(a) (b) 

Figure 3: Examples of diagrams corresponding to the penguin topology: with (a) hard gluon 
and (b) soft gluons. 
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(a) 




Figure 4: Diagrams corresponding to the annihilation topology in the OPE of the correlation 
function KSb^l : (a) factorizable; (b),(c) with hard gluon; (d),(e) with soft gluon. 



Figure 5: Some of the lowest-order diagrams corresponding to the penguin annihilation 
topology. 
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sum rule for the pion decay constant. The diagrams in Fig. [2t»,c describe nonfactorizable 
corrections in the emission topology. The diagrams with soft gluons in Fig. |2fc were cal- 
culated in [13]. Furthermore, penguin contractions for 0\ 2 (some of diagrams are shown 
in Fig. EJ) have also been studied in the framework of LCSR [18] allowing to calculate the 
parameters Tp c , and Tp . In addition, the LCSR for the gluonic penguin operator 

was derived in [17] yielding r^ g n \ In the next two sections, we will present the calculation 
of the remaining diagrams with the annihilation topology (Fig. HJ. 

Having at hand the QCD calculation of a set of diagrams with topology T in terms of 
pion DA's and hard scattering amplitudes, one can then express the correlation function 
in the form of a dispersion relation in the variable s = (p — k) 2 : 

rorw , o 1 7 I™ s Fqcd (s,(p-qf,P 2 ) 

*$8 «p - ^ (p - ^ n = -J *. T-t-ly-i: • (36) 



On the other hand, one can insert a complete set of hadronic states in the 7r meson channel, 
and obtain 



— (p — K) J S — [p — K) 



-, (37) 



where the one-pion ground state contribution contains the pion decay constant and the 
matrix element 



IT£ T ) {{p-q)\P 2 ) =i J d 4 x e-*^<7r-( P -A ; )|r{o(0)j^( a ;)}|7r-(g)> 5 



(3J 



and p\ is the spectral density of heavier hadronic states in this channel. Replacing the 
integral over p\ with the standard duality approximation and equating (|3T)j) to (p?7|). we 
obtain, after the usual Borel transformation: 

((P - d)\ p2 ) = ^~J ds e-^Im s F^ (s + ie, (p - g) 2 , P 2 ) , (39) 

o 

where Sq is the duality threshold in the pion channel. This first step in the derivation 
of LCSR is schematically shown in Fig. |H1 where the diagrams with emission topology are 
chosen for definiteness. 

Note that the hadronic matrix element ()38|) diagrammatically shown in Fig. itself 
represents a correlator. At large negative (p — q) 2 and P 2 it can be factorized into a short- 
distance part (the b-quark propagator) and a long-distance part (the combination of two 
pion DA's). Adding hard-gluon exchanges does not seemingly spoil this factorization. In 
any case, the absence of infrared singularities has to be checked by a direct calculation. 
Hence, instead of using (|3*9*jl. we are free to starting QCD object. That will 
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(a) 



(b) 



Figure 6: The first step in the derivation of the sum rule for B — > tytt amplitude in the 
emission topology: (a) the initial correlation function (only diagrams with emission topology 
are included) is matched to (b) the hadronic dispersion relation in the pion channel where 
only the ground-state pion contribution is shown. 



Figure 7: The second step in the derivation of the sum rule: (a) the two-pion matrix 
element calculated in the spacelike region is analytically continued to (b) the same matrix 
element in the timelike region. 

be done below in the case of the annihilation topology with hard gluons, while employing 
the standard derivation of for the soft-gluon part. 

To continue, following [13] we consider U^' T \(p — g) 2 , P 2 ), as an analytical function of 
the variable P 2 , the invariant mass of the tttc pair. Starting from our calculation for negative 
P 2 , we have to reach the physical timelike point P 2 = m 2 B by analytical continuation (see 
Fig. EJ). As already explained in [13], we use the fact that m 2 B 3> m 2 , Sq, i.e., the two- 
pion system in the B decay is in the timelike asymptotic region, far from the light-quark 
resonances. Note that at fixed (p—q) 2 , the matrix element ni^' 7 ^ can always be represented 
in a form of a hadronic dispersion relation in the variable P 2 : 




(a) 



(b) 



oo 



rt°> T \(p-q)\P 2 ) = lJ 



dt 



Im t n^ J) ((p-g) 2 ,t) 



(40) 



t-P 2 -ie 



o 
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Figure 8: The third step in the derivation of the sum rule: (a) the two-pion matrix element 
H^i' E ^ after analytic continuation to (p — q — k) 2 = m 2 B is matched to (b) the hadronic 
dispersion relation in the B meson channel; the ground state B meson contribution is 
proportional to the B — > tttt matrix element in emission topology. 

This dispersion relation is only needed for illustrative purpose because it allows to unam- 
biguously determine the point of analytical continuation we need as P 2 = m 2 B + ie. Finally, 
it is convenient to represent the QCD calculation result H)^' T \(p — q) 2 ,™^) in the form 
of the dispersion relation in the variable (p — q) 2 , 

nw>«P - ml + fc) = i J* ^r^fc • («) 



and equate this to the hadronic representation in the B meson channel, (see Fig. 8) 



CO , 



where fs is the B- meson decay constant. In the ground-state contribution, we have P 2 = 
(p — k — q) 2 = m 2 B + ie and (p — q) 2 = m 2 B simultaneously, so that the artificial momentum 
k disappears and we encounter the hadronic on-shell matrix element of our interest. After 
applying the duality approximation to the integral over excited states in Eq. (|42|) and 
performing Borel transformation, we obtain the LCSR for the B — » ttti hadronic matrix 
element of a given operator and topology: 

(*-(p)* + (-q)\0\B(p-q)) T = j-^p J ^e^" s W 2 Ims , n (^) (s / + , e)m | + , e) . (43) 
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(a) 



(b) 



Figure 9: Diagrams used to calculate the pion-pion correlator in the annihilation topology. 

Note that the analytical continuations of n^' T ' ) in the variables s' and P 2 interchange 
and may be performed in inverse order as well. 

5 Annihilation with hard gluons 

For the correlation function ()35|) with the operators 0" 2 , the simplest possible diagram 
with the annihilation topology is the factorizable diagram in Fig. |3K- Its contribution is 
expected to vanish due to the conservation of the V — A current for massless u, d quarks. 
In fact, a calculation of the diagram in LCSR yields a parametrically small correction of 
0(sg/m^) which is neglected within the adopted accuracy of the method [13], so that the 
result is consistent with the expectation. 

We start to calculate the annihilation effect by considering the first two 0(a s ) diagrams 
shown in Fig. |3Jd; they contain only the operator 0%. These two-loop diagrams depend on 
four different mass/momentum scales, m&, (p — k) 2 , (p — q) 2 and P 2 , which makes their 
direct calculation technically not feasible. As mentioned in the previous section, instead of 

calculating these diagrams directly, we start from the pion-pion correlator IlS 2 defined 
in Eq. (|38|). with the annihilation topology. The corresponding one- loop diagrams are 
shown in Fig. |U] and we proceed with their calculation. 

(O u A) 

After contracting the quark and gluon fields, the long-distance part of HI,] 2 ' , at leading 
twist 2, reduces to a product of two pion twist-2 DA's: 



where a, (3, 7, 5 and a, 6, c, d are the spinor and color indices, respectively. Convoluting the 
above with the short-distance part, one obtains the following expression for the pion-pion 



(n-( P -k)\d a a (z)u^0) u2(0)d s d (z)\^(l)) 

« (n-(p-k)\d a a (z)u^0)\0) ■ (0\u:(0)d s d (z)\n-(q)) 




1 
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correlator 



{{p-q)\P 2 ) 



where and H% n represent the hard-scattering kernels for the diagrams in Figs. 
and Eb, with the gluon attached to the d and b quark, respectively. The expressions for 
are given in the appendix [B] i n terms of the standard loop integrals. Evaluated at 
(p ~ q) 2 j P 2 < 0) the result for Eq. (|44|) is real and finite, in particular, the numerators in 

Eq. (|44|) cancel at the pole (p — q) 2 u = P 2 v. Furthermore, U^ 2 '^ contains no end-point 
divergences in the variables u and v. The reason why the integrals remain finite is simple: 
the gluon in the diagrams in Fig. El has a virtuality uvP 2 , and remains perturbative unless 
u or v is close to 0, but this region is suppressed by the endpoint behavior of the pion 
DA's. This convergence is important for justifying the replacement of the initial diagrams 
in Fig.|3jD. Note that the latter are convergent and perturbative 'by construction', because 
they contain one of-shell current instead of the pion DA. 

For the sum rule derivation we need the dispersion relation (|4*2*j) for H^ 2 '^ in the 
variable s' = (p — q) 2 . The expression f!44j) continued in this variable, (while keeping P 2 
negative), has branch cuts at s' > 0. The calculation of the imaginary part is straightfor- 
ward, but involved. The result reads: 



a s C F \ f 2 m 



7T 



48 



du / dv 



[(p ~~ l) 2u — P 2 v] 



■ X 



{Ht (% v- (p - q) 2 , P 2 ) + Hi (u, v- (p - q) 2 , P 2 ) } , (44) 



2_2 / U S C F \ f f <p„(u)(p w (v) 



1 1 

du J dv 





\ 1 uv{u s r — v r^y 

s'-^-vP 2 ) 

where u = 1 — u and 



(45) 



v P 2 \m 2 — us'\ 



H x [u, v- s',P 2 ) = uvP 2 (m 2 b -s' + vP 2 ) 2 log (— , , 

\us'{ml — s' + v P Z J 

-(l- ^j) (us' -vP 2 )(m 2 b -us' + vP 2 )(us' + vP 



- (u(m 2 b -us')+v P 2 ) (v m\P 2 + u s'(s' - v P 2 )) log 



s'iuml + uv P 2 ) 



u 



s' 2 + v P 2 (ml — u s') 
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tt I I n 2/ 2 / n2N2, fuvP 2 (ml - s' + V P 2 ) \ 

#2 (u, v; s', P 2 ) = -uv P 2 (m 2 b -s' + v P 2 2 log - - b 2 

\u{s' — vP 2 )\mf —us'\J 

+ (us'-vP 2 ) ( \- , m * \ (m 2 h (us' + (l-2u)vP 2 )-2u(s'-vP 2 )(us'-vP 2 )) 
\ u(s' — V P 2 ) J 

+ (u(m 2 b -us')+v P 2 ) (y P 2 m\ + u s'(s' - v P 2 )) log 



{s' — v P 2 ) (u ml + uv P 2 ^ 



u{u s' 2 + v P 2 {mi — us')) 

The expression (|45jl at negative P 2 contains no singularities within the integration 
region, and is therefore finite. The two pieces proportional to 0(s' — m|) and G(s' — 
m\ju — P 2 v) reflect the two cuts of the diagrams in Fig. El the first one corresponds to 
the on-shell b and d quarks emitted at the B current vertex; the second cut is less trivial 
and emerges when b- and d- quark at the weak vertex are on-shell. 

Finally, according to the procedure explained in the previous section, we analytically 

continue Im s /Ili r < ^ 2 ' A ' ) to the physical timelike point P 2 = m 2 B + ie, so that this function 

(O u A) 

acquires an imaginary part Imp2im s Ti7 r7r 2 ' . The imaginary part in (j45|) naturally origi- 
nates from the logarithms of — P 2 , however the complexity of this expression makes their 
extraction nontrivial. We obtain 



Im p2 Im s ,n^ ( S ',m|) =.^ / > 2 m|K 2 - S ' + ,m|) 2 



x / du j dv -^4^% fe (V -^)-e(s'-m 2 - vm%) 
J J (us' — vm z B y I \ u J v 





(46) 



The existence of nonvanishing imaginary part in P 2 is an important effect we are actually 
looking for, because in the quark-hadron duality approximation it determines the strong 
phase of the hadronic matrix element. Importantly, Eq. (|46J) receives contributions only 
from the diagram in Fig. EK- Physically, the effect corresponds to the d quark from the 
weak decay of b quark going on shell and annihilating with the spectator d quark into a 
virtual timelike gluon. In the diagram in Fig. the gluon is attached to the b quark and 
such mechanism is forbidden kinematically, hence this diagram has no double imaginary 
part. We have verified by explicit calculation that an identical result is obtained if one 
does analytical continuation in P 2 first and only then obtains the dispersion relation in s'. 

Our final result for the annihilation contribution with hard gluons in twist 2 approxi- 
mation is given by the sum rule of the type (|4Hj) : 



(n-n + \O u 2 \B°) h A ard = —^- [ ds' e^- s 'V M ' 2 Im s ,rtSj> A \s',m 2 B + te), (47) 
/ B m B 7r J 

where the real and imaginary part are given by the real part of Eq. ()45|) (with the principal 
value of the integrals containing complex poles) and by Eq. (|46jl. respectively. 
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As already mentioned, the LCSR result (|4T)) is finite, due to the fact that the end-point 
divergence introduced by the gluon propagator 1/uvP 2 is cancelled by the pion DA's. In 
the light-cone expansion of the two-pion diagrams in Fig. |H1 one formally encounters a 
contribution proportional to the two twist-3 pion DA's which is divergent. However, one 
has to keep in mind, that we have used these diagrams only as an effective replacement for 
the part of the diagrams in the correlation function (J33)) . In this function one of the pions is 
interpolated by the axial- vector current, which simply does not have a twist 3 component. 

We conclude that the (twist 3)®(twist 3) contribution to Hl-n- 2 ' has no counterpart in the 
correlation function (j33j) and hence does not play any role in LCSR. 

To complete the calculation of hard-gluon effects, we still need to consider the four 
diagrams in Fig. Ht which belong to the 0(a s ) part of the correlation function (|35|) . Note 
that similar diagrams where B — > tttt annihilation is accompanied by gluons exchanged 
within the weak vertex, also emerge in QCDF approach. However, they were not included 
in [3], where only the diagrams analogous to Fig. were taken into account. The reason, 
apart from expected l/m b suppression, is that in order to describe the two-pion state 
originating from a quark- ant iquark pair, in QCDF one needs an additional long-distance 
object, a sort of two-pion distribution amplitude, which in the local limit reduces to the 
pion form factor at timelike momentum transfer m 2 B . In LCSR approach the 'formfactor- 
like' annihilation diagrams in Fig. emerge as a part of OPE, hence, no new input is 
needed. However, the calculation of these diagrams is not possible with current methods, 
because they contain two loops and many scales. In fact, in this case one cannot use as a 
remedy the pion-pion correlator considered above, since the corresponding diagrams still 
have two loops, even if both pions are described by their DA's. In order to assess this 
effect, remaining at the one- loop level, we employ a completely different method which is 
briefly described in the rest of this section. 

A new type of correlation function is introduced with an on-shell B meson and a pion, 
while interpolating the second pion with a current 5 : 

Il?(p,k,q)=i J d 4 ye^y(0\T{jS(y),d u 2 (0)}\B(p-q)n-(q)). (48) 

We consider the part of this correlation function with quark contractions having annihila- 
tion topology and pick up only diagrams with gluon exchange in the weak vertex. They 
can be obtained from the diagrams in Fig. if the interpolating current is replaced 
by an on-shell 5-meson. The correlation function (|48|) factorizes into a product of vacuum- 
to-pion matrix element (that is, a usual pion DA) and vacuum-to- -B matrix element. The 
latter is expressed via B meson DA's [27,28]. As usual, the external momenta are chosen 
to provide that the virtual quarks and gluon remain far off-shell. For the pion DA we 
retain only twist 2 (twist 3 vanishes in the chiral limit) whereas the two components ± of 
the B meson DA are taken into account. By writing a dispersion relation in the variable 
(p — k) 2 (in the pion channel), one obtains a sum rule for the B — > tttt matrix element. As 

A similar but simpler version of this method was recently applied to B — > 7r form factor in Ref. [29] 
(see also Ref. [30]). 



20 



a result, we find that all four diagrams vanish, which means that contributions from the 
"form factor" part of the annihilation mechanism starts either at higher twists > 4 and/or 
at higher orders in a s and can be neglected. 



6 Annihilation with soft gluons 

The gluons exchanged between the initial 5-meson and the final quark-antiquark state 
in the B — > mi annihilation can also have small virtualities. Diagrams with soft gluon 
contributions cannot be directly calculated in QCDF or PQCD, because in this case it 
is difficult to identify and separate a hard kernel. The soft-gluon nonfactorizable effects 
should either be neglected (arguing that they are l/m& suppressed) or modelled by separate 
nonperturbative parameters. In the LCSR approach the decay amplitude is calculated 
quite differently, by matching the hadronic dispersion relation to the correlation function. 
In the latter, the soft (low virtuality) gluons emerge in OPE diagrams, being emitted 
at short distances and absorbed in the quark- ant iquark-gluon DA's of the pion. For the 
emission topology the corresponding diagrams are shown in Fig. Efc. Their effect, although 
formally 1/rrib suppressed turned out [13] to be of the same order as the 0(a s ) effect 
of nonfactorizable hard gluons calculated from QCDF. On the other hand, the soft-gluon 
effects for the penguin topology (one of diagrams is shown in Fig. Eh) were found suppressed 
in LCSR [18] with respect to the penguin diagrams with hard gluons, indicating that the 
role of soft-gluon effects strongly depends on the topology. Here we will calculate the 
soft-gluon diagrams in the part of the correlation function (|H5jl with annihilation topology. 

The two lowest-order diagrams are shown in Figs. |Ui,e containing an on-shell gluon 
emitted from the heavy-light loop and absorbed in the three-particle pion DA. Technically 
these diagrams are much easier to calculate than the annihilation diagrams with hard 
gluons in FigEJ),c. One returns to the original method of [13] and employs the light-cone 
expansion of the quark propagators in the external gluon field [31]. For the diagram with 
a gluon emission from the massless d quark we use 

i 

S d (x,0) = -i(0\T{d(x)d(0)}\0) = - 2 f 2 . 2 - — ^-r / dvG Tp (vx)(]ta Tp -4ivx Tlp ) + ... , 

2tt 2 (arr \bn z x z I 





whereas the propagator for the massive b quark is 

5^(0,*) = -<<o|r{6(o)5(*)}lo> = / 

i 

" ^ ikx 



(49) 



dvG Tp {vx) I j^-rje , 

(27r) 4 L2 (k 2 — mi) 1 k 2 — m 



m b VXr'Jp 



+ .... (50) 



o 

In the above v = 1 — v, G rp = ^-G^ p and the fixed-point gauge for the gluon field has 
been adopted, having in mind that the hard and soft contributions to OPE are individually 
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gauge-invariant. The dots in Eqs. ()49j) and ()50|) represent terms with derivatives and higher 
orders of the gluon field-strength tensor which we neglect. These terms generate pion DA's 
with twist > 4 and multiplicity > 3, and their contributions to the sum rule are suppressed 
by additional powers of the Borel parameter. 

In the chiral limit for the light quarks, the contribution of the twist-3 quark-antiquark- 
gluon pion DA vanishes and the nonvanishing part comes from the twist 4. The four 
relevant DA's (p\\ ± and (p»± are defined via vacuum-pion matrix elements: 



(0\u(0)t llM G a ^x 3 )d(x 1 )\n(q)) = / w / Va, e -w(*i«i+*»«») 



qz 



(51) 



> 7^75, G af3 - 
Da,i = daida2da^5 (1 — a± — c%2 

where itj are arbitrary numbers and z 



G a/3 and (p^ )± -> y>|| i± . In 
03). The points 2:1,2 are 
is the 



and the one obtained from the above with i^^ 

Eq. (jHU, G a p = \e aPpX GP x 
located on the light cone 
light-cone separation. 

In our case, due to the choice of external spacelike momenta, the regions of integration 
over x, y are close to the light-cone, but x 2 , y 2 , (x — y) 2 are not exactly light-like. Therefore, 
strictly speaking, there is an ambiguity of defining z via x and y in r.h.s. of Eq. (|51J) . From 
the point of view of the light-cone OPE, this ambiguity is a higher twist effect. Indeed, 
calculating the diagrams for different choices z = x, y, x — y in Eq. ()51|) one finds that 
the parts proportional to z a /(qz) in Eq. (JoTj) and its analog for 7^75 yield negligibly small 
contributions. Hence, only two DA's <p± and tp± multiplying the coordinate-independent 
part of the matrix elements appear in the final answer. 

After specifying the propagators and pion DA's, the calculation of the diagrams in 
Figs. IUl,e is straightforward. The following expression is obtained for the invariant ampli- 
tude multiplying (p — k) a in the correlation function : 



(5%, A) _ mff, 

r soft 



da 1 



16vr 2 J -P 2 a x - (p- k) 2 {\- ax) 





(l-ai) 1 1 

da.z I dv 



x 



P 2 (l + 2xv) + 3(p - qf{l - 2vx) 



dx 



m l ~ (p ~ ( ?) 2 (1 — a^v)x 


P 2 (l -Axv) + 3(p-g) 2 



(52) 



Following the procedure explained in section 01 we match the above expression to the 
dispersion relation in the variable (p—k) 2 . After applying duality and Borel transformation, 
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we obtain for the corresponding pion-pion correlator: 



n (02M),soft = j_ m b 



16tt 2 P 2 



111 

ds e~ s/M2 ( da 3 I dv 



dx 



o 











m b ~ (P ~ ?) 2 (1 — Oi 3 v)x 



P 2 {l + 2xv) + 3(p-q) 2 (l-2vx) <p±(0,a 3 ,a 3 )- P 2 (l-4xv)+3(p-q) 2 ^±{0,a 3 ,a 3 ) 



x{l + OK/P 2 )}, (53) 

where a 3 = 1 — a 3 and we neglected terms of 0(s/P 2 ) < 0(sq/P 2 ). 

The analytical continuation to P 2 = m 2 B is then trivial and we can immediately apply 
the dispersion relation in the variable (p — q) 2 , and subsequently, duality and the Borel 
transformation in the B channel. Our final result for the soft-gluon annihilation contribu- 
tion to the hadronic matrix element in the leading twist-4 approximation reads 



it I r>0\soft 



(n-n + \O u 2 \B») 



l6n 2 f B m% 



dse- s ' M2 



ds' 



e m 2 B /M' 2 -s' /M' 2 



1 

du f da 3 



a 3 



u 

2 /.s' 1-u 



m 2 B + 3s' + 2(m 



B 



3s 



V9_l(0, a 3 ,a 3 ) 



m?/s 



ml + Ss'-Am 2 ™ 2 ^-^ 
. B B s'ua 3 



s'ua 3 

0±(O, d 3 , a 3 )] \ x {1 + 0{sl/P 2 )} . (54) 



In the adopted approximation, this part of the decay amplitude does not contribute to the 
strong phase. 

Finally, adding the hard-gluon and soft-gluon contributions given by Eqs. (jj7j) and 
Eq. ([54)1 . respectively, we complete our calculation of the parameter defined in 

Eq. $5). 



7 Factorizable annihilation via Oq operator 

In this section we describe the calculation of the pion scalar form factor defined in Eq. (|29jl . 
According to Eq. ([30)1 . this form factor is needed to estimate the factorizable part of the 
B — > mi matrix element of the operator Oq with annihilation topology. 

In QCDF only part of this contribution was taken into account, namely, the annihilation 
diagrams with the gluon exchange in the final state (see Figs. 4a,b in Ref. [3]). Clearly, 
this is not a complete answer, because the annihilation into two pions via the scalar current 
starts at zeroth order in a s . In addition, there are 0(a s ) diagrams where a hard gluon is 
exchanged in the final state at the vertex of the scalar current. All these contributions lie 
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(a) 



(b) 




Figure 10: Diagrams of the correlation function used to derive LCSR for the pion scalar 
form factor; j = dd is the scalar quark current. 

beyond the usual QCDF approximation: being formally l/m& suppressed, they also do not 
allow a factorization with a hard kernel. One has to parameterize them with a separate 
nonperturbative parameter. 

We use a different approach, considering the scalar pion form factor (n + (p) | qq | n~(q)) 
(q = u, d) as a separate object and obtaining it from LCSR, following the same method 
as in [19]. The calculation is done at spacelike momentum transfer P 2 < 0, analytically 
continuing the result to large timelike P 2 = m 2 B . 

One introduces the correlation function 



where is the same axial- vector current as in Eq. (|H5J) . Inserting the complete set of 
states with the pion and axial meson quantum numbers between the currents in Eq. (|55|h 
one obtains for the ground-state pion contribution: 



At large spacelike p 2 and P 2 = (p — q) 2 the correlation function can be expanded near 
the light-cone and expressed via pion DA's. The OPE starts from the diagram in Fig. ITUk 
which corresponds to the "soft" or end-point mechanism of the pion-to-pion transition. 
In addition, there is a diagram with quark-antiquark gluon DA's (Fig. HUb ) and 0(a s ) 
diagrams, some of them shown in Figs. HOb.d . 




(55) 



T^(p,q) 



^F n s (P 2 ) 

m 2 _ p2 



Pa- 



(56) 
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To calculate the leading-order diagram in Fig. ITUk . the free <i-quark propagator is 
inserted and pion quark-antiquark DA's are factorized out. Only the structure proportional 
to p a is relevant. In the chiral limit (at q 2 = m 2 = 0) the twist 2 contribution vanishes, 
hence the expansion starts from the twist 3. The answer reads: 



T a (p,q) = ifMnf v 



du 



—p 2 u — P 2 u 



(fp{U) 



Id 



Pa 



(57) 



where (pp t(7 are the standard twist-3 DA's of the pion, normalized to fi n = m 2 /(m u +md)- In 
addition, we find that the diagram in Fig. ITHb . with a low-virtuality (soft) gluon absorbed 
in the quark- ant iquark-gluon DA , does not contribute to the relevant Lorentz structure 
~ p a in the chiral limit. 

Furthermore, the experience with the LCSR for the pion vector (e.m.) form factor [19] 
tells that even at rather large momentum transfers of a few GeV 2 the 0(a s ) diagrams 
in Fig. ITUh .d are subdominant. Therefore, in what follows we simplify our calculation 
neglecting their contributions. 

Equating the result ()57|) for the correlation function to the hadronic dispersion relation 
in the variable p 2 , we retain only the pion contribution (JSTIj) and replace the sum over heavier 
states by a quark-hadron duality estimate. For that one has to write the expression ()57|) 
as a dispersion relation in the variable s = p 2 > by substituting the integration variable: 
u = —P 2 /(s — P 2 ). After Borel transformation we obtain the LCSR for the pion scalar 
form factor valid at large spacelike P 2 < 0. Continuing it to large timelike P 2 = m 2 B , we 
get 



[m 



B) 



M ( 1 



B 



1 + OK/mJ) (58) 



u=l 



with an evident end-point dominance. Note that F^(m 2 B ) has no imaginary part in the 
adopted approximation. It will appear at 0(a s ), hence, we expect the strong phase in the 
form factor to be subdominant. 

Substituting Eq. (|HHj) in Eq. (JHUj) . we obtain the estimate for the parameter R^' 6 ' . 
Containing a chiral enhanced factor, and having no a s -suppression, this effect is expected 
to be important. On the other hand, the annihilation effects via 0$$ operators multiply 
small Wilson coefficients C5 and Cq in the amplitude (|27j) . therefore in what follows we only 
retain pS™'^ and neglect nonfactorizable corrections parameterized by r^ r7r ' 6 ' ) and r^ n ' 5 \ 



8 Numerical estimates 

For the numerical evaluation we adopt the same input in all sum rules derived and used in 
this paper, including LCSR (|47|). (}5lj) and (}58|) for the matrix elements in the annihilation 
topology, and the LCSR for the B — >• ir form factor f B7T [16] determining the factorizable 
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amplitude A nn . Finally, fs is substituted by the corresponding two-point sum rule, so that 
the uncertainties of the input parameters partly cancel in the ratios. 

In the LCSR for B — > 7r form factor [16] the one-loop pole mass mj of the b quark is 
used. For consistency, we adopt the same mass for the other sum rules. The current interval 
fhbifhb) = 4.25±0.15 GeV [32] converted into the one-loop pole mass yields mj = 4.7±0.1 
GeV . Furthermore, we take a s (m z ) = 0.1187 [32], evolved to lower scales at two-loop 
order, so that e.g., a s (m b /2) = 0.284. For the factorization and renormalization scale of 
the scale-dependent input parameters, we use //& = a/ m 2 B — m\ ~ 2.4 GeV, numerically 
close to m&/2. Being of order of the Borel parameter M' in the B-meson channel, this scale 
reflects the average virtuality in the correlation function. For consistency, we normalize 
also the Wilson coefficients in H cS at the same scale /j, b . For the latter, we use the NLO 
results obtained in the NDR scheme adopted in our calculation: 



n 


Cl 


c 2 


c 3 


c 4 


c 5 


c 6 


c Sg 


2.4 GeV 


1.124 


-0.272 


0.020 


-0.037 


0.010 


-0.060 


-0.166 


4.8 GeV 


1.073 


-0.174 


0.013 


-0.034 


0.009 


-0.038 


-0.149 



In our uncertainty estimates, we include the variation of the unified scale up to 2/x&. 

Furthermore, for the condensate densities used in the sum rules, we take (<?<?) (1 GeV) = 
(-0.240±0.010 GeV) 3 , so that /^(l GeV) = 1.61±0.20 GeV; (^G 2 ) = 0.012±0.006 GeV 4 
and (qGq) = (0.8 ± 0.2 GeV) x (qq)(l GeV), the scale dependence of the quark-gluon 
condensate density being negligible. The remaining parameters in the B channel are: the 
interval of the Borel parameter M' 2 = 10 ± 2 GeV 2 and the duality threshold Sq = 
35 =f 2 GeV 2 chosen as in [13]. Concerning the parameters related to the pion, we use the 
experimental value / w = 131 MeV [32], and take the Borel interval M 2 = l.Ol^ GeV 2 as 
well as the duality threshold = 0.7 GeV 2 determined from the two-point sum rule for 
U [15]. 

Finally, the pion light-cone DA's deserve a separate discussion. Their definitions, 
asymptotic forms and nonasymptotic parts used in our calculation can be found, e.g., 
in the appendix B of ref. [20] . 

In the twist- 2 pion DA (p n (u) , we include nonasymptotic effects encoded by the Gegen- 
bauer moments a\ and a\. These parameters have recently been estimated in Ref. [35] by 
fitting the LCSR result for B — > n form factor to the data on B — > nlu decay distribution: 

4(1 GeV) = 0.1 ±0.1, 04(1 GeV) > -0.07. (59) 

We combine this range with the constraint [37] obtained from the analysis of the ^7*7 
form factor 

alii GeV) +<(1 GeV) =0.1 ±0.1. (60) 

The resulting intervals 

4(1 GeV) = -h 0.27, aj(l GeV) = -0.07-5- 0.20 , (61) 
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are consistent with other estimates in the literature. Note that Eq. (|5T|) does not exclude 
a purely asymptotic DA. Both Gegenbauer moments are independently varied around 
middle values of the intervals (}6T|) to estimate the corresponding uncertainty. Finally, for 
the twist-3 and 4 pion DA's, the normalization constants and nonasymptotic parameters 
are taken as: / 37r (l GeV) = 0.0035 ± 0.0018 GeV 2 , u 3n (l GeV) = -2.88 and 6*(1 GeV) = 
0.17 ± 0.05 GeV 2 , e w (l GeV) = 0.5, respectively. 

Using the input specified above, we reevaluate the result of [33] for the B —> n form 
factor: 

f+ n (0) = 0.26 ± 0.02 K4] ± 0.03 [pamm] , (62) 

where the uncertainties induced by Gegenbauer moments and by other sum rule parameters 
(the latter are added up in quadrature) are shown separately. This estimate agrees well 
with the most recent calculation of this form factor in ref. [34] where also a small twist-3 
NLO correction is taken into account. In section 2 we added both uncertainties linearly to 
be on the conservative side. In what follows, all errors are added quadratically. 

Using the sum rules for hard and soft gluon contributions presented in sections El and El 
we first estimate the ratio r^ n \ We find that both hard-gluon and soft-gluon annihilation 
contributions, being comparable in magnitude, are numerically very small and partly cancel 
each other, so that 



(-) - [-0.67t{j$ + i (3-6l?;f)] x 10- 3 . (63) 

On the other hand, the factorizable annihilation via quark-penguin operator Oq pro- 
duces a considerably larger hadronic matrix element: 

R { p 6) = 0.231°;^ . (64) 

However, the small Wilson coefficient reduces the effect in the decay amplitudes to the 
same level as for the annihilation via the current- current operator: (c§ + c^/fyR^^ ~ 
Cir^ n \ As mentioned at the end of sect. we therefore neglect both factorizable and 
nonfactorizable 0(a s ) corrections to the annihilation via the operators O5 and Oq. In 
general, the contributions of annihilation amplitudes are found at the same level as the 
other nonfactorizable effects estimated from LCSR in Refs. [13,17,18]. 

Having at hand the new estimates of and R^' 6 ^, we now update the phenomeno- 
logical analysis of B — > tttc channels, with all nonfactorizable parts of the amplitudes cal- 
culated from LCSR, except the emission with hard gluons which is estimated using QCDF. 
To this end, we recalculated the previous LCSR predictions using the current input, which 
has only slightly changed. For the nonfactorizable emission, we obtain 



r E 



(7T7r) 



(1.818;|) xKT 2 + (1.3113 xlO- 2 + z(-4.7l^)xKT 2 , (65) 



- soft 



- hard 



where the soft-gluon part is obtained from [13] and the hard-gluon contribution is estimated 
using QCDF [2,3] with the default value and error for the parameter describing the twist-3 
hard-spectator diagrams. For the quark-penguin operators with scalar-pseudoscalar Dirac 
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a$ r P {B + -> 7T+7T ) 


a%p(B° -> 7T + 7T-) 


ag£(fl° -> 7T 7T ) 


BaBar 


-0.01 ±0.10 


0.09 ±0.16 


0.12 ±0.56 


Belle 


0.02 ±0.08 


0.56 ±0.14 


0.44 ±0.56 


Average 


0.01 ±0.06 


0.37 ±0.10 


0.28 ± 0.40 


This work 





-0.04 ±0.01 ±0.01 


n 7n +0. 19+0.08 
u - ' u -0. 29-0.08 



Table 1: Direct CP-asymmetries from experiment (all numbers taken from [25]; the errors 
added in quadrature) compared with the LCSR predictions. 



structure, we have found that the soft-gluon emission contribution vanishes in twist < 4, 
and for the hard-gluon part we again use QCDF: 

r^ 6) = [(-2.7±0.4)xl0- 2 ] hard . (66) 

The penguin-topology effects are calculated from LCSR obtained in Refs. [17,18]. The 
resulting ratios to the factorizable amplitude are: 

r { ;; ] = [o.ni£g ± i (i.i±S:?)] x io- 2 , r ™ = [-0.1818S + < (-o.soi® ] x io- 2 , 

r£° = (0.93±g:g|) x IO" 2 , r£*> = - (3.8^) x 10~ 2 , (67) 

and the modified penguin parameters are fp^ ~ rp*' — ( tt sCV)/(367r); for brevity we do 
not show the corresponding numbers. 

Having specified all parameters entering the decomposition of B — > tttt amplitudes in 
Eqs. ()2())1 and ()21)1 . we calculate the branching ratios and direct CP-asymmetries, using 
the values of B meson lifetimes from [25], and the relevant CKM parameters from [22]. In 
particular, we adopt \V ub \ = (4.22 ± 0.26) ■ 10~ 3 (the errors added in quadrature) and use 
a representative interval 7 = (58.6 ± 10)°. The results are: 

BR{B + - tt+tt ) = (6.7±i;|±g;l) xlO" 6 
BR(B° - 7T+7T-) = (9.7±?;i±i;l) xlO" 6 

BR(B°^n\°) = (0.291°;^;°?) xl(T 6 , (68) 

where the errors represent the variation of the LCSR parameters and of the CKM factors, 
respectively. The direct CP asymmetries are presented in Table ^ For completeness, we 
also calculate the amplitudes T, P and C defined in f)34jl parameterizing them as in [5]: 

«' A = ^, de^ = -^. (69) 

-*- TT7T -*- 7T7T 



We obtain: 



x = 0.29±g;5| , A = (-21^)° 



d = 0.22^ : ™ , = (-173±l) , (70) 
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where the second error in d stems from the uncertainty of \V ub \. 

We find that the general picture does not qualitatively deviate from the naive factoriza- 
tion considered in sect 2. Although the nonfactorizable contributions substantially enhance 
the small B° — > 7r°7r° amplitude predicted in naive factorization, the disagreement between 
theory and experiment for the branching ratio of this decay remains. Our prediction (J70J) 
for the amplitudes is inconsistent with the fit [5] to data. As already discussed in sect. 2, 
one possible interpretation of this disagreement is a missing contribution to the 1 = parts 
of the amplitudes T nn , C nn and P n7r . 



9 Comparison with QCD factorization 

Let us first investigate the behavior of the annihilation amplitudes obtained from LCSR 
in the heavy quark limit, by making standard substitutions in the sum rules: 

m B = m b + A, Sq = m 2 b + 2m b u , M' 2 = 2m b r, f B = fs/y/rr^, (71) 

where A,u>q,t and fs are m b -independent parameters 6 . At m b — * oo, the factorizable 
amplitude A n7T defined in Eq. (I12|) scales as m l J 2 . Expanding in l/m b the annihilation 
contribution with hard gluons given by Eq. (|37j) and dividing it by A^, we reduce the real 
part of the diagram in Fig. (gluon emitted by the light quark) to 

i i 

rA ^~m b j dU —j dV — + - (?2) 



where A is a generic energy scale not related to m b and higher powers in l/m b are denoted 
by ellipses. Although formally suppressed by l/m b , this expression recovers the logarithmic 
divergence at u — > 0, of the annihilation contribution in QCDF. Importantly, in LCSR this 
divergence occurs only at m b — > oo. Moreover, as the numerical analysis in sect. |H1 shows, at 
finite m b ~ 5 GeV the logarithms originating from the end-point region (at m b — > oo they 
reduce to Log(m b /r) and hog{m b /uj)) do not produce an enhancement of the annihilation 
contribution. Furthermore, analyzing LCSR in the heavy quark limit, we observe that the 
phase generated by the diagram in FigEK, as well as the contribution of the gluon emission 
from the heavy quark (FigEb), remain finite at m b — > oo, being of the same 0(l/m b ). 
Other contributions that we have calculated from LCSR, the soft-gluon part of and 
the factorizable annihilation with Oq given by R^ n ' 6 \ are also finite, being suppressed by 
an additional power of l/m b with respect to Eq. (|72|l . 

The origin of the end-point divergence in QCDF and the reason why it is absent in the 
LCSR approach can be readily understood. The power counting in QCDF implies that the 
momentum of the light quark in the B meson has to vanish, i.e. the propagator of the light 

6 Note that the scale r has to be large in comparison with Aqc d , allowing one to use the power 
expansion in Aqcd/t, so that QCD sum rules remain valid in the heavy quark limit (the well known 
example is the sum rule for fs in HQET [36]). 
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quark shrinks to a point [2,3]. This is different in LCSR, since the B meson is effectively 
replaced by the spectral density of the heavy-light quark loop integrated over the duality 
interval ml < s' < Sq , and thus the momentum of the light quark is non- vanishing. In 
other words, the endpoint singularity is regularized by the typical momentum of the light 
quark. 

To demonstrate this, let us slightly modify QCDF by convoluting the annihilation hard 
kernels with the standard .B-meson DA. We consider the annihilation diagrams in FigEl 
where instead of the current an on-shell B meson is inserted, represented by its 

DA. The result for the annihilation amplitude A\ (see the definition in [3]) in the twist 2 
approximation reads: 



A\ = 7ra s J du(p^(u) J dmp^iu) J dvip T (v) 
ooo 



x 



1 u + u/rriB 



uv(u — u/rriB^ uv(l — (u — uu/m B )(v — u/m B )j 



(73) 



where the .B-meson DA 4> B (u) is normalized as: J °° ducj)~^(uj) = 1. Neglecting the spectator 
quark momentum u, that is, replacing (f) B (oo) — > S(uj), one recovers the expression for A\ 
given in [2,3] (see also [38]), with an end-point divergence in the first term in brackets 
(corresponding to the diagram with gluon emission from the light quark). However, it is 
generally expected (see e.g., [28]) that (fi B (uj) — > at uj — > 0. The integral in Eq. (|73*|) 
then converges (taken as a principal value) yielding Log(m B / \ B ), \ B being the size of the 
region in u where the function <P b {uj) dominates. Simultaneously, this expression acquires 
an imaginary part, 7 due to the pole in the integration region at u = uo/m B . 
Employing a realistic model of (p B (uj), e.g., the one suggested in [28]: 

= ^~ wIXb , (74) 

it is easy to calculate Eq. (|73*|) numerically and to estimate the corresponding parameter 
r^ r7r ' ) (see Appendix A for the relation between this parameter and A\). The result turns 
out small, with both real and imaginary parts at the level of ~ 1%, that is, roughly of the 
same size as the LCSR estimate obtained in the previous section. The model of annihilation 
represented by Eq. (fT3*|) is rather crude, because the transverse momenta of the quarks in 
B meson are neglected, but their account could not qualitatively change the result. Thus, 
from the point of view of the LCSR approach, the end-point divergence in the annihilation 
diagrams in QCDF originates in the hard-scattering approximation and the m& — > oo limit. 
In the phenomenological analysis of B — > tttc done in Ref. [3], a model for the annihi- 

7 The appearance of an imaginary part due to the momentum of spectator quark was noticed already 
in [6] while discussing the differences between QCDF and PQCD. 
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lation diagrams was used, replacing all divergent integrals by a generic logarithm: 

1 

[^^X A =(l + p A e^)\n^, (75) 

J y A h 



where = 0.5 GeV, pA < 1 and the phase (p A is arbitrary. With this model the effective 
annihilation coefficients {B^/A^bi, (where B wn = i(Gp/^/2)f B f%) entering the decay 
amplitudes have been estimated [3]. Using our results and the relations given in Appendix 
1X1 we can also obtain these coefficients: 

[-0.15^ + ^(0.82^)] xlO- 2 , b 2 = % u 

[-1.31& + i (0.015±g5S)] x 10- 2 , h - . (76) 

ci 

The above estimates, especially for b\, differ from the ones presented in [3] for the default 
value pa = in Eq. (|75|). In fact, a numerical agreement is not anticipated, because the two 
sets of fej's originate from two different methods. Moreover, one cannot expect that LCSR 
predictions for fej allow a parameterization with a single complex parameter X A - Important 
is that both the default values of the annihilation coefficients in [3] and our predictions in 
Eq. (|7f)jl are very small in comparison with the factorizable amplitude. Hence, LCSR is in 
a qualitative agreement with QCDF, if the annihilation effects in the latter are represented 
by moderate logarithms. 

To complete the comparison with other methods, let us briefly discuss PQCD [6]. In 
this approach, all B — > amplitudes as well as the B — > 7r form factor are represented by 
the diagrams with 0(a s ) hard-scattering kernels and meson wave functions. Another dis- 
tinctive feature of PQCD concerns nonvanishing transverse momenta of partons in mesons, 
which can only be introduced in a model-dependent way. Hence, the annihilation ampli- 
tudes in this approach are protected from the end-point divergences and acquire imagi- 
nary parts. There is however a basic difference between LCSR and PQCD approaches to 
B — > ixix. In LCSR the diagrams at 0(a s ) are subleading and numerically suppressed, 
justifying the perturbative expansion within OPE. Moreover, the higher-twist soft-gluon 
diagrams of 0(1/ m^) are as important as the 0(a s ) effects. Importantly, the dominant 
part of the B — > 7r form factor in LCSR is "soft" and has no relation to a s . In PQCD, the 
whole form factor (hence, the factorizable part of the B — > tttt amplitude) is assumed to 
be of 0(a s ). Furthermore, the main contribution to the strong phase in PQCD stems from 
the annihilation mechanism with the scalar-pseudoscalar operator Oq, and the diagrams 
again start at 0(a s ) level. We have also found Oq to be an important source of factoriz- 
able annihilation, but in LCSR the 0(a s ) contribution to this mechanism is expected to 
be subleading in comparison with the zeroth order in a s , "soft" contribution which has 
been calculated in sect.0 Starting at 0(a s ) level and neglecting soft contributions, PQCD 
nevertheless predicts annihilation effects that are larger then in LCSR. This surprising fact 
means that it is difficult if not impossible to reconcile these two approaches with each 
other. 



B m 



A 

B, 



7T7T 
7V7C ■ 



A 



31 



10 Conclusion 



In this paper, the weak annihilation contributions to B — > ttty decay amplitudes have 
been calculated applying the method of LCSR. This work complements previous stud- 
ies [13, 17, 18, 40] of nonfactorizable effects in B — > nn with the QCD sum rule approach. 
In LCSR, due to sufficient virtuality of the underlying correlation function, the OPE dia- 
grams with annihilation topology are free from end-point divergences. Both contributions 
of hard and soft gluons are taken into account. A finite result for the hadronic matrix 
element of the current-current 0\ operator with annihilation topology is obtained includ- 
ing an imaginary part which contributes to the strong phase. In addition, an important 
factorizable contribution from the quark-penguin operator 6 has been found. For the 
annihilation with hard gluons considered in this paper, we have modified the method sug- 
gested in [13], to avoid the problem of calculating two-loop multiscale diagrams. Instead 
of performing the QCD calculation based on the vacuum-to-pion correlation function, we 
start from the pion-pion correlator, thereby reducing the calculation to one-loop diagrams. 

We emphasize that QCD sum rules have a limited accuracy, at the same time one is able 
to estimate the uncertainties of the method. Moreover, many uncertainties cancel in the 
ratios of nonfactorizable and factorizable hadronic matrix elements obtained from LCSR. 
Obtaining B — > tttt hadronic matrix elements one uses an additional assumption of the 
local quark-hadron duality, allowing the transition (analytical continuation) from a large 
spacelike scale to the large timelike scale m B . This approach has much in common with 
evaluating the timelike asymptotics of the pion e.m. form factor from the QCD calculation 
in the spacelike region. 

Our main phenomenological result is a smallness of the annihilation contributions in 
B — > 7T7T. This is generally consistent with QCDF, if the divergent annihilation diagrams 
there are modelled by moderate logarithmic factors. Hence, we find no compelling reason 
to consider the annihilation amplitude as a free parameter in QCDF. The relatively large 
values of this parameter generated by the fits to the data are probably not originating from 
the annihilation mechanism. 

Small annihilation effects predicted in this paper are in the same ballpark as other 
nonfactorizable contributions obtained from LCSR, including charming and gluonic pen- 
guin topologies and nonfactorizable corrections to the emission topology. Altogether, the 
smallness of the corrections to the leading-order factorizable amplitude reveals a good con- 
vergence of the OPE series for the correlation function and justifies the use of the adopted 
approximation in LCSR, that is, including only 0(a s ) and twists < 4, as well as omitting 
the small 0(sl/m B ) corrections in each term of OPE. 

Furthermore, we have performed the phenomenological analysis of the three B — > iyty 
channels using the results of LCSR for all nonfactorizable effects, except the hard-gluon 
nonfactorizable corrections to the emission topology. For the latter the default prediction 
of QCDF [2,3] is used. Our results disagree with the current data for BR(B° — > 7r + 7r~) 
and BR(B° — > -r°-k°) and probably also for the direct CP-asymmetry in B° — > ir + ir~, that 
is, for the channels where the AI = 1/2 weak transition (or, equivalently, the 1 = two- 
pion final state) contributes. If the experimental data do not change, we have to admit 
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that LCSR misses an important part of the isoscalar amplitude. Then a new "inverse 
AI = 1/2" rule has to be established for these decays, meaning that the amplitude A Q 
introduced in sect. 2 has to decrease after including the missing piece. 

A natural question arises: are there additional mechanisms in B — > tttt which may fill 
this gap. The penguin-annihilation topology comes first to one's mind, a possibility which 
was not yet explored by both QCDF and LCSR. This effect contains multiloop diagrams 
and cannot be easily evaluated. Still, the experience with LCSR for charming penguins [18] 
tells us that the OPE diagrams of the type shown in Fig. El could not be large. In any 
case, the penguin-annihilation mechanism deserves a closer look in future. One can for 
example speculate about nonperturbative gluonic effects in this mechanism, which have an 
anomalously large scale (3> Aq C d) and lie beyond OPE. Such effects will most probably 
influence other neutral final states, e.g., enhance B — > p°p°. That however seems not to 
be the case because the current experimental upper bound for BR(B — > p°p°) is less than 
the measured value of BR(B — > 7r°7r°). 

Another resource of enhancement is the hard-spectator part in the nonfactorizable 
emission. The corresponding diagrams in QCDF diverge at twist 3 and are replaced by 
another generic logarithm. In our numerical estimates we have taken the default value 
of this parameter from [3], but the fits to the current data [4] with a free parameter for 
this contribution produce large effects. Indications that the hard-spectator mechanism 
is important, were found recently in the SCET framework [39]. We plan to study the 
nonfactorizable hard-gluon emission, including the hard-spectator mechanism in LCSR 
(work in progress). 

Finally, an important avenue of future studies is the sum rule analysis of the B — > ttK 
and B — > KK channels including calculable SU(3) violation effects which are expected [40] 
to be important. The variety of kaon channels with a lot of accumulated data will allow 
to isolate various topologies, and to put the LCSR approach under a tighter scrutiny. 
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A Effective coefficients in QCDF 

For convenience, we present here the relations between the nonfactorizable matrix elements 
parameterized in section El by the ratios and the effective coefficients dj introduced 
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in [2,3] for B — > tttt. The decay amplitudes written in terms of these coefficients are: 
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where r£ = 2fi n /mb an d electroweak penguins are neglected. Comparing with Eqs.(|20|). 
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The correspondence between the annihilation diagrams introduced in [3] and parameters 
is schematically given by 
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where the sign = indicates that in LCSR and QCDF different approximations are used to 
calculate l.h.s and r.h.s., respectively. In particular, R^' 5 ^ is of the zeroth order in a s , 
whereas A{ is of 0(a s ). 
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B Kernels for hard annihilation 

The kernels for the hard annihilation contribution to the correlation function in Eq. f)44j) 
are given by 

(u,v;s' = (p-q) 2 , P 2 ) = \p 2 (2m 2 b - s'(2 + u) + 3P 2 v) 1 B (P 2 uv, 0, 0) 
- 1 

-UV 
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B and Co are the standard two-point and three-point functions, respectively: 
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